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We consider a 5-dimensional Ricci flat bouncing cosmological model in which the 4- 
dimensional induced matter contains two components at late times - the cold dark matter 

■ (CDM)+baryons and dark energy. We find that the arbitrary function f(z) contained in 
the solution plays a similar role as the potential V(<j>) in quintessence and phantom dark 
energy models. To resolve the coincidence problem, it is generally believed that there is 
a scaling stage in the evolution of the universe. We analyze the condition for this stage 
and show that a hyperbolic form of the function f(z) can work well in this property. We 

Q ' find that during the scaling stage (before z fa 2), the dark energy behaves like (but not 

| identical to) a cold dark matter with an adiabatic sound speed c 2 s fa and p x 0. After 

^ , z S3 2, the pressure of dark energy becomes negative. The transition from deceleration 

■ to acceleration happens at zj< fa 0.8 which, as well as other predictions of the 5D model, 
agree with current observations. 
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1. Introduction 



In recent decades, the observations of high redshift Type la supernovae have re- 
vealed that our universe is undergoing an accelerated expansion rather than decel- 
erated expansion 1 - 2 < 3 . Meanwhile, the discovery of Cosmic Microwave Background 
(CMB) anisotropy on degree scales together with the galaxy redshift surveys in- 
dicate fltotai — 1 4 and Q m ~ 1/3. All these results strongly suggest that the 
universe is permeated smoothly by dark energy' which has a negative pressure and 
violates the strong energy condition. The dark energy and accelerating universe 
have been discussed extensively from different points of view 5 > 6 ' 7 . In principle, a 
natural explanation to the cosmic acceleration is the cosmological constant. How- 
ever, there exist serious theoretical problems such as the fine tuning problem and 
coincidence problem. To overcome the fine tuning problem, a self-interacting scaler 
field, dubbed quintessence, with an equation of state (EOS) = p<f> / Pc/> was in- 
troduced, where is time varying and negative. Thus, by properly choosing the 
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forms of the potential V (0), desired behaviors of the quintessence can be obtained: 
(i) with a negative pressure which drives the universe accelerated expansion; (ii) 
with an energy density that was much smaller than that of the ordinary matter 
(and radiation) at early times (due to the constrains from primordial nucleosyn- 
thesis and structure formation) and is comparable to the latter at recent times. To 
resolve the coincidence problem, tracker solutions or scaling solutions were designed 
in which there exists a stage where the ratio of potential and kinetic energies of the 
scalar field remains constant approximately. In this way, all the initial 'information' 
could be eliminated. So, in the past, if the ratio of the densities of dark energy and 
matter is approximately constant, the coincidence problem could be avoided. 

It has been drawn great attention to the idea that our universe is a 4-dimensional 
hypersurfce embedded in a higher dimensional world as is in Kaluza-Klein theories 
and in brane world scenarios. Here we consider the 5-dimcnsional Space- Time- 
Matter (STM) theory 11 . This theory is distinguished from the classical Kaluza- 
Klein theory by that it has an non-compact fifth dimension and that it is empty 
viewed from 5D and sourceful viewed from AD. That is, in STM theory, the 5D 
manifold is Ricci-flat which implies that the five-dimensional space-time is empty, 
the matter of our universe is induced from the fifth dimension, and the AD hy- 
persurface is curved by this induced matter. Mathematically, this approach is sup- 
ported by Campbell's theorem 12 which says that any analytical solution of Einstein 
field equation of N dimensions can be locally embedded in a Ricci-flat manifold of 
(N + 1) dimensions, though there is an argument recently concerning Campbell's 
theorem 13 . 

In this paper we consider a class of five-dimensional cosmological models 9 of 
the STM theory. This class of exact solutions satisfies the 5D Ricci-flat equations 
Rab = and is algebraically rich because it contains two arbitrary functions fi (t) 
and v (t). It was shown 9 that several properties characterize these 5D models: (i) 
The AD induced matter could be described by a perfect fluid plus a variable cosmo- 
logical 'constant', (ii) By properly choosing the two arbitrary functions, both the 
radiation-dominated and matter-dominated standard FRW models could be recov- 
ered, (iii) The big bang singularity of the AD standard cosmology is replaced by a 
big bounce at which the universe reaches to a finite and minimal size. Before the 
bounce, the universe contracts; after the bounce, it expands. Also, the evolution 
of the universe containing three components was discussed in 10 . In this paper we 
assume that the universe is permeated smoothly by CDM+baryons p m as well as 
dark energy p x with pressure p x = w x p x at late time. We will show that a scaling 
solution can be obtained which is different from the conventional scaling solutions 
and late-time attractors: (i) There exists a scaling stage before the transition from 
decelerated to accelerated expansion, and so the coincidence problem could be re- 
solved; (ii) During this stage, the dark energy has the equation of state w x ~ and 
behaves like a cold dark matter - this would provide a unified model for dark matter 
and dark energy as the Chaplygin gas model does 15 . Given an appropriate form of 
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f(z), we will find that before z w 2, the dark energy behaves like a 'tracker' solution 
and mimics CDM+baryons. The transition from deceleration to acceleration is at 
zt « 0.8, which agrees with recent cosmological observations. 

This paper is organized as follows. In section 2, we assume that the induced 
matter of the universe contains two components: CDM+baryons and dark energy. 
We write down Einstein field equations and the equation of states (EOS) of dark 
energy. In section 3, we rewrite these equations via redshift z, then the arbitrary 
functions /i(t) contained in the 5D solution corresponds to a function / (z), which 
is found to play a similar role as the potential V{4>) in quintessence and phantom 
models. With use of a kind of hyperbolic function of / (z), a scaling stage is obtained. 
Also, we discuss the evolution of the dimensionless density parameters, the EOS of 
dark energy, the deceleration parameter and the clustering of dark energy. Section 
4 is a conclusion. 



2. Dark energy in the 5D model 

Within the framework of STM theory, a class of exact 51? cosmological solution 
was presented by Liu and Mashhoon in 1995 14 . Then, in 2001, Liu and Wesson 
9 restudied the solution and showed that it describes a cosmological model with a 
big bounce as opposed to a big bang. The 5D metric of this solution reads 

ds 2 = B 2 dt 2 - a 2 (y^2 + rW ) - d v 2 C 1 ) 



where dQ 2 = (d0 2 + sin 2 6d<f> 2 ) and 

A 2 = (n 2 + k)y 2 + 2uy 



v 2 + K 
V 2 + k' 



Here fi = /z(f) and v = u(t) are two arbitrary functions of t, k is the 3D curva- 
ture index (k = ±1, 0), and K is a constant. This solution satisfies the 5D vacuum 
equation Rab = 0. So, the three invariants are 

h = R = 0,I 2 = R ab Rab = 0, 

h = RabcdR ABCD = ^ 8 ■ (3) 

The invariant -Z3 in Eq. (3) shows that K determines the curvature of the 5D 
manifold. Here we should mention that because the 5D field equations of the STM 
theory are always sourceless, i.e., Gab = or Rab = 0, there is no need to introduce 
a higher-dimensional Newtonian constant or a higher-dimensional Planck mass in 
the theory. This is one of the main differences between STM theory and other 
Kaluza-Klein theories. 
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Using the 4D part of the 5D metric (1) to calculate the 4D Einstein tensor, one 
obtains 

( 4)r o _ 3 (M 2 + k) 

G - ^2 ' 

2/i/t p 2 + k 

In the previous work 9 , the induced matter was set to be a conventional matter plus 
a time variable cosmological 'constant'. In this paper, we assume that the induced 
matter contains two parts: one is CDM+baryons p m with pressure p m = 0, the 
other is the dark energy p x with pressure p x . So, we have 

3 (p 2 + k) 



WG\ = =< 4 > Gl = ^ + (4) 



A 2 



= Pm + Px, 



AA + A 2 - Px > (5) 

where 

Px = w x p x . (6) 
From Eqs.(5) and (6), one obtains the EOS of dark energy 

Px 2pp/AA + (p 2 + k)/A 2 
Wx Px 3(v 2 + k)/A 2 -p m0 A-3> [t) 
and the dimcnsionless density parameters 

Pm+Px~ 3(p 2 +k)A> ^ 
Q x 1 ^m- (9) 

where p m o — p m o^o- The proper definitions of the Hubble and deceleration param- 
eters should be given as 9 , 

, . J 2 A I fdA\ 2 Ap 

from which we see that fi / p > represents an accelerating universe, (i / p < 
represents a decelerating universe. So the function p{t) plays a crucial role in deter- 
mining the properties of the universe at late times. It was pointed out that another 
arbitrary function v (t) may relate closely to the early epoch of the universe 10 . 



3. Late time evolution of the cosmological parameters versus 
redshift 



We consider the spatial flat case k = 0. From the solution (2) we see that on a 
given y = constant hypersurface we have A = A(t). So, without loss of generality, 
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one can always write \i — u(z) and v = v{z) if y is fixed, where z is the redshift, 
A I A = 1 + z. Now we define u 2 E j y? z = f (z) and then we find that Eqs. (7)-(ll) 
reduce to 

l + (l + z)dlnf(z)/dz 
3-3n m0 (l + z)f(z)/f(0)' [ ' 

f(z)(l + z) _ (l + z)f(z) 

1 + Z E J (V) 

n =i-n + (w) 

t + z E 

_ 1 + m x w x (l + z)d\nf{z) 
Q - 2 _ ~~ 2 di - " (15) 

Here the subscript E denotes the point where fl mZE — Q XZE (which is 1/2 when the 
radiation is neglected). As is known from quintessence and phantom dynamical dark 
energy models, the potential V ((f) is unfixed. By choose different forms of V (</>), 
desired properties of dark energy could be obtained 17 . Now, there is an arbitrary 
function / (z) in the present 5D model. Different choice of / (z) corresponds to 
different choice of the potential V [4>). This enables us to look for desired properties 
of the universe via Eqs. (12)-(15). 

To obtain a scaling solution, there should be a stage where the ratio of the 
energy densities of dark energy and dark matter is approximately constant. So, the 
condition for exiting a scaling stage at late time is jf^ ~ const. From Eqs. (13) and 
(14), this condition leads to 

^<«H 1 + «>Uq,. (16) 
dz 

Also, to obtain an accelerating universe at late times, some mechanics is needed 
to break the scaling state steeply and make the dark energy dominated. In 18 , a 
hyperbolic function k((p) = k m i n + tanh(4> — tfii) + 1 as a leap kinetic term to achieve 
the desired features of quintessence is discussed. So, in our case, to have a scaling 
dark energy, we let f(z) being of a similar form, 

/w = tanh[ f + ; b)]+c , d7) 

where, a, b and c are arbitrary constants and can be estimated by observation data, 
i.e. by using Q. mZB = 0.5, n m0 = 0.273±^, q Q = -0.67 ± 0.25. Then we obtain 
the values of ze, a and b as 

1 (1 + 6) (1 + A) z E ln±±£ 
a -2^ n (l-<5)(l-A)' b -^ (i + *)(i + A) - ( 18 ) 

where, 5 = 1 + z E — c and A = c — (1 + ze) Q m0 ■ Using the equation qo = 
— i d ln d ^ z ^ , we can solve ze numerically by leaving c as a free parameter which 
could be determined by the constraints from structure formation (In this paper we 
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-1.41 



take c = 1.45). So, the function / (z) describes the evolution of CDM+baryons and 
dark energy completely at late time. 

Under the choice of the function f (z) in the form of (17), the ratio of the 
dimensionless density parameters f2 m and Q x becomes constant when z > z s , that 
is to say, the dark energy is scaling with the matter. And the EOS of dark energy is 
w x w 0, though it is negative in nature. Thus we conclude that just not a long time 
ago, the EOS became negative and the dark energy began to drive the universe 
to accelerate. For different values of tt m o an d qo, the parameters Zt, ze and w x o 
are listed in Table 1. We can see that the transition point from deceleration to 
acceleration is at zt ~ 0.8, which is compatible with current observations. Also, 
for a larger f2 m o, a larger absolute value of w x q is needed to drive the universe 
accelerating. For the case of fl m o = 0.232, q = —0.45 (which corresponds to 
a = 4.52 and b = 0.29), we plot the evolution of dimensionless density parameters 
fl m , £l x , EOS of dark energy w x , and deceleration parameter q versus redshift z in 
Fig. 1. The adiabatic sound speed of the dark energy is 

2 _ i '-'/'• / ( ... dPx dw x 



c ' = w *-te V°*-te- p '-te*- (19) 



where, in terms of the redshift z, the pressure p x is 

/4(l + z) 2 f (l + z)dlnf(z) 



Alf(z) V dz 

-oV(o)^(<™a +1 ). 

We plot the adiabatic sound speed (? s in Fig. 2, where a = 4.52, b = 0.29 and 
c = 1.45. From Fig. 2 we find that the adiabatic sound speed of the dark energy 
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Fig. 1. The dimensionless density parameters versus redshift z: Q x of dark energy (dotted line, 
the third line from the top), Q m of CDM+baryons (solid line, the first line from the top), EOS 
of dark energy w x (dashed line, the forth line from the top), and the deceleration parameter q 
(dash-dotted line, the second line from the top). The figure is plotted with a = 4.52, b = 0.29 
and c = 1.45. In these values, the current values of the cosmological parameters are f2 m o = 0.232, 
&x0 = 0.768, w x o = —0.82 and qo = —0.45, respectively. From z = 10 to z = 400, we find 
/ttm ~ 0.02 which is indeed a scaling. 

is c 2 rts during the scaling stage. So, it behaves like the conventional matter with 
an almost zero w x . However, for z < 2, the sound speed c 2 s becomes wild and even 
negative - this is also encountered in other unified dark matter and dark energy 
models. Note that before z — 2 the dark energy with c 2 s s=s and w x w behaves like 
a dark matter (but not identical to it because w x is negative in nature). However, 
one should take into account the clustering of dark energy seriously. As did in Ref. 
20 , we consider a universe containing only the dark energy with c 2 s w and the 
perturbation 5 = Sp x /p x , where p x is the dark energy density. Since w x ^ 0, the 
Jean's instability equation is modified, 

d 2 5 3 , r 9 . d8 ( k 2 c 2 3 , 9 ,\ „ , . 

where a is the FRW scale factor, H = — (the prime denotes the deriva- 
tive w.r.t. the conformal time), ^[c 2 ,^] = 1 — 5w x + 2c 2 and B[c 2 ,w x ) = 
1 — 6c 2 + 8w x — 3w 2 . In general, the "growing" solution is S oc a 7 where 7 = 
i (1 - \A + [(1 - § A) 2 + 6B] 1 / 2 ). For c 2 = w x = 0, 7 = 1, the solution S oc a is 
the conventional dust-dominated growing solution. For a small but none-zero c 2 
and w x , the dark energy may cluster during scaling. However, the process is slower 
than the conventional dust-dominated case. And once w x < —0.12 (for all c 2 > 0), 
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Fig. 2. The evolution of the adiabatic sound speed versus the redshift z. The figure is plotted 
with a = 4.52, b = 0.29 and c = f.45. In these values, the current values of the cosmological 
parameters are fi m o = 0.232, Q x o = 0.768, w x o = —0.82 and go = —0.45, respectively. 

7 will less than zero and the dark energy stop clustering. As 19 does, the growth 
exponent / is defined as 

_ d In 4(a) . , 

d In a 

which is roughly k- independent for a wide range of k. And, the growth factor g of 
density perturbations between arbitrary a\ < a 2 is defined as 

5 k (a 2 ) ( a 2 \ f 

5(fll ' a2) "4M = Uj (23) 

with a suitably defined average growth exponent /. For example, the scaling stage 
begins at the dark matter dominated point adec (or Zd ec = 1100, decoupled epoch) 
and ends at a es (or z = 2), we have the growth factor g(adec, a es ) w 1100. After the 
scaling stage and if there is no clustering as expected, we have g(a es , a ) ~ z( s and 
/ = — log Zes (zdec — z es ) = — log 2 (1097) = —10 roughly. In our case, the evolution 
of 7 versus redshift z is plotted in Fig. 3 numerically. In Fig. 3, we find that the 
singular points exist which correspond to a negative infinity for 7. So, after this 
stage, we cannot find any clustering dark energy in the universe at present. However, 
we must point out that the wild properties of the dark energy might be out of our 
knowledge and make the dark energy more mysterious. 

The influence of early dark energy on CMB anisotropy was discussed in 21 . It 
was pointed out that large Sl x is ruled out by cosmic observations. In our case, 
the dimensionless energy density Q x in scaling stage depends on the parameter 
c mainly, we can choose the parameters properly to fit observational constraints. 
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Fig. 3. The evolution of 7 versus redshift z. In the scaling stage, / ~ 1 behaves like dark matter. 
In the breaking stage, three singular points exist. The figure is plotted with a = 4.52, b = 0.29 
and c = 1.45. In these values, the current values of the cosmological parameters are U m o = 0.232, 
f^xO = 0.768, WxO = —0.82 and qo = —0.45, respectively. 



From Fig. 1, we find that Cl x w 0.02 at scaling stage which is in compatible with 
observations. In 19 , a general formula to calculate cr 8 was given as 



~ (£*!)3(n'V)-n;' (s))/5 (1 _ noyJz 1 -^ 1 )^ M A ) ^ (24) 
a s (B) a eq ' Y t (B) 

In our case, during the scaling stage, the dark energy behaves like dark matter 
with a zero pressure. So, comparing the model with the standard cold dark matter 
(SCDM) (ct s {SCDM) rj 0.5 ± 0.1), we have 

<j s (x) = 1100- 3(0 - 02 - 0)/ V 8 (S*(7L>M) = 0.37cr 8 (SCDM) w 0.46. (25) 

Moreover, observations such as cluster abundance constraints yield 22 

<j$ = (0.5±0.1)r! m 7 (26) 

where 7 is slightly model dependent and usually 7 w 0.5. So, in our case, we have 

af us -{x) = 0.5 x 0.98~ - 5 sa 0.51. (27) 

So, the ratio is 

- 0-92, (28) 
which is close to unity and compatible with observations. 
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4. Conclusion 

A general class of 5D cosmological models is characterized by a big bounce as 
opposed to the big bang in 4D standard cosmological model. This exact solution 
contains two arbitrary functions fj,(t) and u(t), which are in analogy with the dif- 
ferent forms of the potential V ((f)) in quintessence or phantom dark energy models. 
Also, once the forms of the arbitrary functions are specified, the universe evolution 
will be determined. In this bounce model, we assume the universe to contain at 
late time two components: CDM+baryons and dark energy. Instead of choosing the 
forms of the arbitrary functions \x (t) and v (t) , we transform \i (t) from coordinate 
time t to redshift z. So, the choice of /z (i) becomes the choice of / (z). Then, we 
let / (z) to be of the form of a hyperbolic function (17) and study the evolution of 
the universe. We find that there is a scaling stage before z « 2, in which the dark 
energy has constant ratio with CDM+baryons and behaves like a pressureless cold 
dark matter with w x « 0. After that, the pressure of the dark energy deviates from 
w x w and gradually becomes much more negative. The transition from deceler- 
ated expansion to accelerated expansion is at Zt ~ 0.8 which, together with other 
predictions of the present 51? model, is in agreement with current observations. 
Thus we see that the arbitrary function / (z) of the 5D model plays a similar role 
as the potential V (</>) in quintessence or phantom models in eliminating the initial 
'information', then the coincidence problem is resolved. Also, this kind of model 
would provide a unified model for dark matter and dark energy. 

However, we also see from Fig 2 and 3 that several singular points for z < 2. 
Perhaps, these singularities correspond to some unknown phase transition. 

Acknowledgments 

This work was supported by NSF (10273004), (10573003) and NBRP 
(2003CB716300) of P. R. China. 

References 

1. A.G. Riess, et.al., Observational evidence from supernovae for an accelerating universe 
and a cosmological constant, 1998 Astron. J. 116 1009, astro-ph/9805201; S. Perlmut- 
ter, et.al., Measurements of omega and lambda from 42 high-redshift supernovae, 1999 
Astrophys. J. 517 565, astro-ph/9812133. 

2. J.L. Tonry, et.al., Cosmological Results from High-z Supernovae , 2003 Astrophys. J. 
594 1, astro-ph/0305008; R.A. Knop, et.al., New Constraints on CIm> ^A, an d w from 
an Independent Set of Eleven High-Redshift Supernovae Observed with HST, astro- 
ph/0309368; B.J. Barris, et.al., 23 High Redshift Supernovae from the If A Deep Sur- 
vey: Doubling the SN Sample at z^O.l, 2004 Astrophys. J. 602 571, astro-ph/0310843. 

3. A.G. Riess, et.al., Type la Supernova Discoveries at z > 1 From the Hubble Space 
Telescope: Evidence for Past Deceleration and Constraints on Dark Energy Evolution, 
astro-ph/0402512. 

4. P. de Bernardis, et.al., A Flat Universe from High- Resolution Maps of the Cosmic Mi- 
crowave Background Radiation, 2000 Nature 404 955, astro-ph/0004404; S. Hanany, 



February 5, 2008 6:48 WSPC/INSTRUCTION FILE scaling'manu 



Scaling Dark Energy in a Five- Dimensional Bouncing Cosmological Model 11 

et.al., MAXIMA-1: A Measurement of the Cosmic Microwave Background Anisotropy 
on angular scales of 10 arcminutes to 5 degrees,2000 Astrophys. J. 545 L5, astro- 
ph/0005123; D.N. Spergel et.al., First Year Wilkinson Microwave Anisotropy Probe 
(WMAP) Observations: Determination of Cosmological Parameters, 2003 Astrophys. 
J. Supp. 148 175, astro-ph/0302209. 

5. I. Zlatev, L. Wang, and P.J. Steinhardt , Quintessence, Cosmic Coincidence, and the 
Cosmological Constant, 1999 Phys. Rev. Lett. 82 896, astro-ph/9807002; P.J. Stein- 
hardt, L. Wang , I. Zlatev, Cosmological Tracking Solutions, 1999 Phys. Rev. D 59 
123504, astro-ph/9812313; M.S. Turner , Making Sense Of The New Cosmology, 2002 
Int. J. Mod. Phys. A 17S1 180, astro-ph/0202008; V. Sahni , The Cosmological Con- 
stant Problem and Quintessence, 2002, Class. Quant. Grav. 19 3435, astro-ph/0202076. 

6. R.R. Caldwell, M. Kamionkowski, N.N. Weinberg, Phantom Energy: Dark Energy 
with w < —1 Causes a Cosmic Doomsday, 2003 Phys. Rev. Lett. 91 071301, astro- 
ph/0302506; R.R. Caldwell , A Phantom Menace? Cosmological consequences of a 
dark energy component with super-negative equation of state, 2002 Phys. Lett. B 545 
23, astro-ph/9908168; P. Singh, M. Sami, N. Dadhich, Cosmological dynamics of a 
phantom field, 2003 Phys. Rev. D 68 023522, hep-th/0305110; J.G. Hao, X.Z. Li , 
Attractor Solution of Phantom Field, 2003 Phys. Rev. D 67 107303, gr-qc/0302100. 

7. Armendariz-Picon, T. Damour, V. Mukhanov, k-Inflation, 1999 Physics Letters B 458 
209; M. Malquarti, E.J. Copeland , AR. Liddle, M. Trodden, A new view of k- essence, 
2003 Phys. Rev. D 67 123503; T. Chiba , Tracking k-essence, 2002 Phys. Rev. D 66 
063514, astro-ph/0206298. 

9. H.Y. Liu and P.S. Wesson, Universe models with a variable cosmological "constant" 
and a "big bounce", 2001 Astrophys. J. 562 1, gr-qc/0107093; T. Liko, P.S. Wesson, 
gr-qc/0310067; S.S. Seahra, P.S. Wesson, Universes encircling five-dimensional black 
holes, 2003 J. Math. Phys., 44 5664; Ponce de Leon J, 1988 Gen. Relativ. Gravit. 
20 539; L.X. Xu , H.Y. Liu, B.L. Wang, Big Bounce singularity of a simple five- 
dimensional cosmological model, 2003 Chin. Phys. Lett. 20 995, gr-qc/0304049; H.Y. 
Liu, Exact global solutions of brane universe and big bounce, 2003 Phys. Lett. B 560 
149, hep-th/0206198; B.L. Wang, H.Y. Liu , L.X. Xu , Accelerating Universe in a Big 
Bounce Model, 2004 Mod. Phys. Lett. A 19 449(2004), gr-qc/0304093. 

10. L.X. Xu and H.Y. Liu, Three Components Evolution in a Simple Big Bounce Cosmo- 
logical Model, astro-ph/0412241, to be published in IJMPD. 

11. P.S. Wesson Space-Time-Matter (Singapore: World Scientific) 1999; J.M. Overduin 
and P.S. Wesson, Phys. Rept. 283, 303 (1997), gr-qc/9805018. 

12. J.E. Campbell, A Course of Differential Geometry, (Clarendon Oxford, 1926); S. 
Rippl, R. Romero, R. Tavakol, 1995 Gen. Quantum Grav. 12 2411; C. Romero, R. 
Tavako and R. Zalaletdinov, 1996 Ge. Relativ. Gravit. 28 365; J. ELidsey, C. Romero, 
R. Tavakol and S. Rippl, 1997 Class. Quantum Grav. 14 865; S.S. Seahra and P.S. 
Wesson, Application of the Campbell-Magaard theorem to higher-dimensional physics, 
2003 Class. Quant. Grav. 20 1321, gr-qc/0302015. 

13. E. Anderson, The Campbell-Magaard Theorem is inadequate and inappropriate as 
a protective theorem for relativistic field equations, gr-qc/0409122; F. Dahia and C. 
Romero, Dynamically generated embeddings of spacetime, gr-qc0503103. 

14. H.Y. Liu and B. Mashhoon, A machian interpretation of the cosmological constant, 
1995 Ann. Phys. 4 565. 

15. F. Pedro Gonzalez-Diaz, Unified Model for Dark Energy, astro-ph/0212414; V.F. Car- 
done, A. Troisi, S. Capozziello, Unified dark energy models : a phenomenological ap- 
proach, Phys.Rev. D 69 (2004) 083517, astro-ph/0402228; M.C. Bento, O. Bertolami, 
A. A. Sen, The Revival of the Unified Dark Energy-Dark Matter Model ?, Phys.Rev. D 



February 5, 2008 6:48 WSPC/INSTRUCTION FILE scaling'manu 



12 Lixin Xu, Hongya Liu 

70 (2004) 083519, astro-ph/0407239; Hongsu Kim, Brans-Dicke Theory as an Unified 
Model for Dark Matter - Dark Energy, astro-ph/0408577; Zong-Hong Zhu, Generalized 
Chaplygin gas as a unified scenario of dark matter/ energy: observational constraints, 
Astron.Astrophys. 423 (2004) 421, astro-ph/0411039. 

16. H.Y. Liu, Exact global solutions of brane universe and big bounce, 2003 Phys. Lett. B 
560 149, hep-th/0206198. 

17. V. Sahni, Theoretical models of dark energy, 2003 Chaos. Soli. Frac. 16 527. 

18. A. Hebecker and C. Wetterich, Natural Quintessence ?, 2001 Phys. Lett. B 497 281. 

19. M. Doran, J.M. Schwindt, and C. Wetterich , Structure formation and the time de- 
pendent of the quintessence, 2001 Phys. Rev. D 64 123520. 

20. S. DeDeo, R.R. Caldwell, P.J. Steinhardt, Effects of the Sound Speed of the 
Quintessence on the Microwave Background and Large Scale Structure, Phys. Rev. 
D 67 (2003) 103509, astro-ph/0301284; T. Padmanabhan, Structure Formation in the 
Universe, Cambridge University Press (1995). 

21. R.R. Caldwell, M. Doran, CM. Mueller, G. Schaefer, C. Wetterich, Early 
Quintessence in Light of WMAP, 2003 Astrophys.J. 591 L75. 

22. L. Wang and P.J. Steinhardt, Astrophys. J. 508, 483 (1998). 



